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Abstract. The multiwavelength emission spectrum and 
associated velocity field of the Seyfert prototype NGC 
4151 is modeled. NGC 4151 has been since the thirties sub- 
ject of extensive modeling, nuclear photoionization being 
the basic approach considered by all authors. HST data 
has impressively revealed the existence of a large range 
of velocities (100 - 1500 kms~^) dominating the emitting 
clouds in the extended emission line region of the galaxy. 
Following this observational result, a revision of the pho- 
toionization modeling approach applied to NGG 4151 is 
presented. It is concluded that a mixture of radiation dom- 
inated clouds and shock dominated clouds are required to 
explain the multiwavelength line and continuum spectra 
of the galaxy. The relative contribution of shock excitation 
versus photoionization is consistently modeled along the 
nebulae taking into account the spatial variation of both 
flux and velocity of main optical and UV lines. The mul- 
tiwavelength continuum spectrum of NGG 4151 is then 
nicely accounted for by a combination of nuclear emission 
at high energy, gas Bremsstrahlung, and dust emission. 
The last two phenomena are directly linked to the compos- 
ite shock and photoionization excitation of the gas. The 
radio SED is found however dominated by synchrotron 
emission created by Fermi mechanism at the shock front. 
In addition, a 3x10'^ K black body component accounts for 
the host galaxy contribution. As a result of the modeling, 
silicon is found depleted by a high factor and included in 
dust grains, while N/G abundance ratio is found compat- 
ible with cosmic values. 

Key words: galaxies : nuclei - galaxies : Seyfert - shock 
waves - galaxies : individual : NGG 4151 - X-ray : galaxies 



1. Introduction 

Detailed imaging of the narrow-line region (NLR) of 
Seyfert galaxies with the Hubble Space Telescope (HST) 
has revealed its complex morphology and velocity field. 
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Recently, observations of NGG 4151 indicated the pres- 
ence of emitting clouds with velocities ranging from -1-846 
to -1716 kms"i(Winge et al. 1997, Kaiser et al. 1999). 

NGG 4151 is a nearby barred galaxy, usually classified 
as a Seyfert type 1, although it has already been con- 
sidered as a Seyfert 1.5 (Osterbrock & Koski 1976) and 
has also shown the characteristics of a Seyfert 2 galaxy 
(Pension & Perez 1984). NGC 4151 is one of the most ob- 
served active galactic nuclei (AGN), from radio to X-rays. 
The recent HST longslit emission-line data (Nelson et al. 
2000), coupled to ISO data (Alexander et al. 1999), offer 
an excelent opportunity to improve our understanding of 
the NLR physical conditions with the added effect of the 
velocity field. 

The kinematics derived from the long slit observations 
(Nelson et al. 2000, Crenshaw et al. 2000) show evidence of 
three components: a low velocity system, consistent with 
normal disk rotation, a high velocity system in radial out- 
fiow at a few hundred kms~^and an additional high ve- 
locity system with velocities up to 1400 kms~^, as pre- 
viously found from STIS slitless spectroscopy (Hutchings 
et al. 1998, 1999, Kaiser et al. 1999). The authors see 
the signature of a radial outflow, with no interaction with 
the radio jet. However, high velocity components, shifted 
up to about 1500 kms~^ from the systemic velocity, are 
also seen by Winge et al (1997) associated with individual 
clouds located preferentially along the edges of the radio 
knots. Such association suggests a cloud-jet interaction, 
which deflnitively may influence the morphology and the 
physical conditions of the NLR. 

The high spatial resolution permits to see that the 
emission line ratios vary substantially on scales of a few 
tenth of an arcsecond, indicating that the density and ion- 
ization state of the emitting gas are strongly influenced by 
the local conditions, hence suggesting that shock fronts 
may be at work. 

The presence of high velocity clouds in the NLR of 
active galaxies has been predicted from the modeling of 
the emission-line and continuum spectra of the Seyfert 2 
objects NGC 5252 and Circinus (Contini, Prieto & Viegas 
1998a, 1998b), as well as suggested from the more gen- 
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eral analysis of the optical-ultraviolet continuum of a rela- 
tively large sample of Seyfert 2 galaxies (Contini & Viegas 
2000). High velocity clouds are also invoked in an alter- 
native model based on a mnltiwavclcngth analysis of the 
emission-line and continuum spectra of the source (Con- 
tini & Viegas 1999), which explains the soft X-ray excess 
in NGC 4051. 

In all of these models the physical conditions of the 
high velocity clouds are usually shock-dominated. The 
shocks contribute to the high ionization emission- lines and 
to the continuum emission, mainly by providing efficient 
heating of dust, contributing to reproduce the observed 
middle-infrared emission, and to the observed soft X-ray 
spectrum which originates in the high temperature post- 
shock zone. 

Yet, the interpretation of the NGC 4151 emission-line 
and continuum spectra is so far residing exclusively on 
photoionization models. All the authors agree that, given 
the simplicity of the models, the observed line ratios can 
be reproduced. 

In NGC 4151 the high velocity clouds are directly ob- 
served. Velocities of ^ 1500 kms^^ could be theoretically 
explained extrapolating to the NLR the arguments about 
high velocity clouds raised by Weymann et al. (1997). 
Here, however, in order to obtain a large range of velocities 
collisions are invoked. 

Although photoionization models reproduce within 
better than a factor of 2 the most important line ratios 
(Schulz 1988, 1990; Kraemer et al. 2000; Alexander et al. 
1999), some questions remain still open, as for example 
the mere fact of the influence of the velocity field in the 
cloud spectrum. Thus, in this paper we will look for the 
signature of the high velocity gas in the emission-line spec- 
trum, as well as in the radio to X-rays continuum spectral 
energy distribution (SED). 

We present a new type of single-cloud modeling based 
on the spatial distribution of the observational data in 
the optical range; then, a multi-cloud model is proposed 
to explain the emission lines in the IR. 

As was done previously in modeling other objects 
(NGC 5252, Circinus), the multi-cloud model issued from 
the analysis of the emission lines is then constrained by fit- 
ting the continuum SED in a large frequency range. This 
is possible only with composite models (shock -|- photoion- 
ization), because the gas heated only by the radiation from 
the active center cannot reach temperatures high enough 
to fit the Bremsstrahlung emission and dust reradiation 
in the large range of the observed frequencies. 

A review focused in the successes and problems of pho- 
toionization models is presented in Sect. 2. The compos- 
ite models, coupling photoionization and shocks, are de- 
scribed in Sect. 3. The optical emission-line spectrum is 
discussed in Sect. 4, while the infrared lines are presented 
in Sect. 5, and the results for the ultraviolet lines appear 
in Sect. 6. The observed and calculated SED are compared 
in Sect. 7, and the conclusions appear in Sect. 8. 



2. A brief review of photoionization models 

Most recently, photoionization models for the NLR of 
NGC 4151 have been proposed by Alexander et al. (1999) 
and Kraemer et al. (2000). In both cases, the observed 
emission-line ratios are reproduced within a factor of 2. 
Both set of authors use, however, different modeling ap- 
proach: the matter distribution adopted is different. In 
Alexander's et al. (1999) models (hereinfter called type 
A), the best fit is found with a single cloud component 
with a filling factor less than unity; in Kraemer's et al. 
(2000) models (type B), a multi- cloud is used, namely a 
less dense matter-bound component is added to a dense 
radiation-bounded component. We note that Alexander 
et al. discarded this latter solution arguing that a matter- 
bounded component implies three free additional param- 
eters, whereas the introduction of a filling factor -their 
case- implies only one. Nevertheless, let us recall that 
a matter-bounded component favors the high ionization 
lines, while (radiation-bounded) models with filling-factor 
less than unity mimic a lower average density, leading to 
a larger ionized zone and favoring the low-ionization lines. 
Thus, models assuming a filling factor less than unity are 
not a good physical representation of a clumpy zone; yet, 
the spatial distribution of [SII] ratios across NGC 4151 ex- 
tended emission line region reveals instead a rather clumpy 
region with regions of higher and lower density at different 
points along the HST slit (Nelson et al. 2000). Regarding 
both type A and type B models, the authors agree that 
"the fit of the line ratios is good taking into account the 
simplicity of the models" . 

Simplicity has always been a strong argument in fa- 
vor of photoionization models applied to nebular regions 
in AGNs, in addition to the undisputable presence of a 
strong central radiation source. However, in Science it is 
usually from the attempt to explain the "imperfections" 
of a model, (those data not explained by it), that a more 
realistic scenario can be drawn. With that in mind we list 
below the discrepancies between the above proposed pho- 
toionization models and the observational data of NGC 
4151. 

Due to the various difficulties to have a self- consistent 

data set (data from different epoch, different resolutions 
and apertures), Alexander et al. (1999) use several criteria 
to select the set of emission-lines to be reproduced by type 
A models. One of the criteria excludes lines emitted by 
ions that "can be easily ionized by other processes". The 
set of emission-lines excluded includes some usually used 
in the diagnostic diagrams of nebular gas: He II 4686, [O 
III] 4363, [O I] 6300, [N II] 6548-^6584 and the density 
indicator [SII] 6716, 6731 doublet, in addition to the high- 
ionization Fe lines: [Fc X] 6734 and [Fe XI] 7892. The final 
set of lines used for modeling, that includes UV, optical, 
and IR lines, is reproduced within a factor of 2 by their 
best fit models. Discrepancies larger than two are derived 
for [Fe VII] 5721, [Fe VII] 6086, [Ne II] 12.8 and [Ne III] 
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36.0 which arc undcrprcdicted, and for [S III] 9069, [Ne 
V] 14.3, [S III] 18.7, [Ne V] 24.3, and [S III] 33.5 that are 
overpredicted. 

Regarding the SED imphed by models, the presence of 
a big blue bump in the ionizing radiation spectrum, peak- 
ing at 50 eV, is excluded by their best-fit model because 
its effect would be to overproduce the low- ionization lines 
and underproduce the high-ionization lines. However, this 
is the net effect produced by imposing a filling factor less 
than unity, as adopted in their best-fit model. 

Kraemer et al. models (type B ) use HST/STIS low- 
dispersion long-slit data at position angle P.A. = 221°. 
Thus, their data set includes UV and optical lines, but 
no infrared lines that are the basis of type A models by 
Alexander et al. (1999). Their model results are compared 
with observed emission-line ratios derived at different po- 
sitions along the slit. The large majority of the lines are 
reproduced within a factor of 2. Divergencies larger than 2 
are often found for the high ionization lines: [Fe VII] lines, 
and the UV lines N V 1240, [Ne IV] 2423, and NIV] 1486, 
although in the latter case the signal-to-noise is very low. 
At some locations of the emitting nebulae, the calculated 
C III] 1909, CII] 2326, [O III] 4363, and the [S III] lines 
also largely diverge from the data. Notice, however, that 
CIII] and CII] are blended and that, in the case of the [S 
III] lines the discrepancy may be due to an instrumental 
effect as pointed out by Kraemer et al. (2000). 

Several interpretations for the above discrepancies are 
discussed by Kraemer et al. (2000). Our assumption is 
that those discrepancies may be revealing the presence of 
an additional ionizing mechanism, which is not the domi- 
nant process, but that shows its signature through partic- 
ular observational features. This point of view is adopted 
here assuming that the additional mechanism is due to the 
presence of shocks. 

3. Single-cloud models 

Faint high velocity emission regions intermingled with 

brighter emission clouds are shown in NGC 4151 imag- 
ing (Hutchings et al. 1999). We consider this observational 
fact as an indication that the extended NLR of NGC 4151 
is a mixture of low velocity radiation-dominated clouds 
and high velocity shock-dominated clouds, all contribut- 
ing to the emission-line spectrum. Accordingly, composite 
models accounting for the coupled effect of the central ion- 
izing radiation and shock excitation due to cloud motions 
are assumed. Numerical simulations for single clouds are 
obtained with the SUMA code (see, for instance, Viegas & 
Contini 1994). Notice that the simulations apply whether 
the shocks originate from an interaction of the emitting 
clouds with the radio jet or from a radial outflow of the 
clouds. 

The input parameters are the shock velocity, Vs, the 

preshock density, no, the prcshock magnetic field. Bo, the 
ionizing radiation spectrum, the chemical abundances, the 



dust-to-gas ratio by number, d/g, and the geometrical 
thickness of the clouds, D. A power-law, characterized by 
the power index a and the flux, Fh, at the Lyman limit, 
reaching the cloud (in units of cm^^ s^^ eV^^) is gener- 
ally adopted. For all the models, Bq = lO"** gauss, auv 
= 1.5, and ax =0.4, and cosmic abundances (Allen 1973) 
arc adopted. The basic models are calculated with d/g = 
10~^^, however, this value is changed a posteriori to better 
flt the continuum SED. 

Shock dominated models (SD) are calculated assuming 
that the effects of the shock prevail on radiation (Fh=0.). 
Radiation dominated models (RD), however, are compos- 
ite, i.e. they account both for photoionization and shocks 
up to Vs=500 kms^^, but photoionization dominates the 
physical conditions of the emitting gas. 

The grid of models which arc actually used for mod- 
eling is presented in an accompanying paper by Contini 
& Viegas (2001, hereafter referred to as CVOl). In the 
following, calculated emission-line ratios from a selected 
number of models in the grid are compared to the HST 
log-slit optical data at P.A.= 221° (Nelson et al. 2000 and 
Kraemer et al. 2000) and to the ISO integrated aperture 
SWS data by Sturm et al. (1999). Models are selected on 
the basis of the physical conditions of the emitting gas 
dictated by the observations, e.g. the FWHM of the line 
profiles for Vg. 

In models which account for shock effects, the density 
downstream is determined by compression, which depends 
on the shock velocity, and changes considerably with dis- 
tance from the shock front (CVOl, Figs. 5a, 6a, and 7a). 
Therefore, the preshock density and the shock velocity 
chosen define a distribution of the density across the cloud, 
which must be adequate to provide a good fit to the den- 
sity sensitive lines, e.g., [Oil] 3727, [NI] 5200, [Nil] 6548 
([SII] 6716-1- is not very significant because S can be locked 
in dust grains). 

Moreover, in the NLR the density of the clouds fol- 
lows the gradient of the cloud velocity, decreasing with 
the distance to the center. For each cloud, the observed [S 
II] 6717/6730 line ratio (which does not depend on S/H) 
is used as a first test for the choice of no and Vg. Then, 
the intensity of the ionizing radiation, the physical condi- 
tions calculated by the model, and the geometrical thick- 
ness of the cloud are deduced from the line spectrum, as a 
whole. The [O III]/H/3 line ratio is an indicator for Fh. On 
the other hand, wide clouds are optically thick, leading to 
stronger low ionization level lines, while narrow clouds are 
matter-bounded, with fainter low-ionization lines. Thus 
the choice of D is then constrained by the best fit of a 
large number of line ratios, particularly, the ratio of the 
low ionization lines to H/3. 
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4. The optical emission-lines 

Modeling the spatial distribution 
Because of the availability of spatially resolved emission- 
line spectra in different regions of NGC 4151, a new mod- 
eling approach is followed in this work. The most signif- 
icant optical emission lines are modeled accounting for 
their spatial distribution across the NLR. 

In each position wc check the consistency of the pre- 
vailing models explaining several line ratios, including 
those that are poorly reproduced by the photoionization 
models. 

The observed and calculated emission-line intensities 
relative to H/3 are presented in a series of figures (Figs. 1- 
4), where the emission- line ratios are shown as a function 
of the projected nuclear distance, including both the SW 
data (on the left) and the NE data (on the right) at P.A.= 
221°. 

In order to determine the intensity of the power-law 
radiation flux in the NLR edge closer to the nucleus a 

preliminary estimate is made by comparing model results 
with the observed line ratios as presented in Fig. 1. Filled 
squares refer to Nelson et al. 2000 (not reddening cor- 
rected) and filled triangles to Kraemer et al. 2000 (red- 
dening corrected) data, respectively. 

The scales for log Fjj are indicated upon the upper hor- 
izontal axis. Each curve represents RD models calculated 
with different Fh but corresponding to the same Vg- ng : 
Vs=100 kms-i and no=100 cm'^ (dotted lines), Vs=200 
kms^i and no=2GG cm~^ (dash-dotted lines), Vs=300 
kms~^ and no=3GGcm~^ (short-dash lines), Vs=500 and 
no=300 cm~^, (long-dash lines). Thin lines refer to a nar- 
row cloud (D=1G^''~^^ cm) and thick lines to wider clouds 
(D = 1019 cm). 

A sequence of three different cases (log Fh=13, 12, 
and 11.3) is shown in the top, middle, and bottom di- 
agrams, respectively. We have chosen the [OIIIJ/H/3 and 
[0II]/H/3 line ratios, because are the most significant. The 
fit of both [0III]/H/3 and [0II]/II/3 data is acceptable only 
in the bottom diagrams (notice that in these diagrams the 
flux is not fully symmetric, indicating that the SW and NE 
regions are slightly different). So, for consistency we chose 
the bottom case to model also the other line ratios. 

Since we are investigating the intensity of ionizing ra- 
diation from the active center (AC) and the velocity dis- 
tribution in the central region of NGC 4151, both Fh and 
Vg are shown in Figs. 2-4 diagrams. In this way, the rela- 
tive importance of the velocity field and the AC radiation 
field can be recognized. The scales for log Fh and Vs(in 
kms^i) are indicated upon the horizontal upper axis of 
the top and middle and bottom diagrams, respectively. 

In the top diagrams of Figs. 2-4 (panels a and d) each 
line corresponds to RD model results as in Fig. 1. Panels 
b and e show the solid lines corresponding to SD mod- 
els (CVOl, Tables 1-10) for which the maximum value of 
the velocity distribution is Vs=7G0 kms"^ in the central 



region, while diagrams c and f correspond to the results 
with a maximum velocity of Vs=14GG kms"^. 

Because shock dominated models for narrow and wide 
clouds give very similar results, these models are repre- 
sented by one line (model results overlap), corresponding 
to one serie of results. We draw attention to the thin and 
thick lines in the two bottom diagrams which refer to re- 
duced and full intensity line ratios, respectively, and not 
to models calculated by small and large D, as shown in 
the top diagrams. 

The model results shown in Figs. 2-4 cover most 
of the data consistently for the [OIII] 5007+4959/H/3, 
[Oil] 3727/H/3, [01] 6300+6363/H/3, [SIII] 9G69/H/3, [SII] 
6713-F6731/H/3, and [FeVII] 6086/H/? hue ratios. The first 
three ratios are related to the ionization state of the gas, 
while the other three are chosen because they are not well 
reproduced by pure photoionization models. 

The ranges and the spatial distributions of Fh were 
chosen phenomenologically by the fit, particularly, of the 
[OIII]/H/? and [SIII]/H/3 ratios, while the ranges and the 
distributions of Vg were chosen by fitting the low ioniza- 
tion {[Oil] /up, [SII]/H/?) and neutral ([0I]/H/3) line ra- 
tios. The data corresponding to the high ionization line 
ratio [FeVII] /H/3 are not enough to constrain the models. 

We adopt a velocity field decreasing from the center of 
the galaxy outwards. The observed velocity field is, how- 
ever, complex, so, two different shock velocity distribu- 
tions arc shown, one with a maximum 700 kms"'^ in 
the central region (middle diagrams), and another with 
a higher maximum Vg, up to 1400 kms"^ (bottom dia- 
grams) . 

Two different lines represent the SD models in both 
middle and bottom diagrams. The thick one corresponds 
to the calculated models, the thin one to the models down- 
wards shifted by a factor of 3 (middle diagrams) and of 
2 (bottom diagrams). The shift in the middle diagram is 
dictated by the fit of the [0I]/II/3 ratios (sec Sect 4.1) 
and in the bottom diagram by the [FeVII] /H/3 ratios (see 
Sect. 4.3). These shifts do not represent lower abundances 
of the elements relative to H, but they indicate that in a 
multi-cloud model corresponding to the weighted sum of 
single-cloud models, the weights of the SD high velocity 
(> 300 kms"'^) clouds is reduced (see Sect. 5). Indeed, by 
modeling the data on a large scale, this solution may look 
arbitrary, particularly considering faint lines (e.g. [01]), 
and lines affected by the presence of dust (e.g. [SIII], [SII], 
etc.). If reducing the weights of the SD models, the fit 
of all the emission-lines is consistently improved, we can 
conclude that this reduction is soimd and SD models cor- 
responding to high velocities in the nuclear region have 
lower weights. The relative weight accounts for the rel- 
ative number of clouds in the conditions determined by 
the model, for the dilution factor (i.e., the square ratio 
between the distance of the cloud from the galaxy cen- 
ter) and the distance of the galaxy from Earth. Notice 
that the middle and bottom diagrams in Figs. 2-4 refer 
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Fig. 1. The [O III] 5007-h4959/H/3 (left diagrams) and the [O II] 3727/H/3 (right diagrams) emission-line ratios are 
shown as a function of the projected nuclear distance. 




Fig. 2. The [O III] 5007+4959/H/3 (left) and the [O II] 3727/H/3 (right) emission-line ratios are shown as a function 
of the projected nuclear distance. Same notation as in Figure 1. 




Fig. 3. The [O I] 6300+6360/H/3 (left) and the [S III] 9069/H/3 (right) emission-Hne ratios as a function of the projected 
nuclear distance. Same notation as in Figure 1. 




Fig. 4. The [S II] 6713+6731/H/3(left) and the [Fe VII]6086/H/3 (right) emission-line ratios as a function of the 
projected nuclear distance. Same notation as in Figure 1. 
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to "high velocities" in the nuclear region of the galaxy 
(between 1500-700 kms~^ancl about 300 kms~^). Indeed, 
Nelson et al (2000, Fig. 5) show that the bulk of cloud 
velocities is within 300 kms^'^. So, the reduction of the 
weights in the SD diagrams indicates that the number of 
high velecity SD clouds is small. 

Recall that observational data result from the integra- 
tion along the line of sight, which may include different 
clouds. Thus, a more realistic fit should be obtained by a 
multi-cloud average (see Sect. 5). A compromise bcitwecn 
a consistent picture over a large spatial distribution of 
several emission lines and the precision of the fit must be 
achieved. 

Observed features 

In several cases there arc two or more emission lines 
within the passband, so that different velocities are sam- 
pled for the different lines (Hutchings et al. 1999). These 
are the [Oil] 3727 and [SlI] doublets, as well as [Nil] lines 
plus Ha at about 6548 A. The [OIII] image covers veloc- 
ities between -1200 and -860 kms~^. High velocity mate- 
rial in this velocity range is much fainter than the main 
bright clouds showing low velocities. It is seen on both 
sides of the nucleus and outside the main biconical emis- 
sion regions. The high velocity gas is weak in H/3 and [OH] 
lines and seems to be associated with highly ionized ma- 
terial. Ionization of oxygen is higher along radial locations 
on both sides of the nucleus. Generally, there is an asso- 
ciation of high velocity clouds with high ionization gas, 
however, there are regions of high ionization (mainly to 
the E side) with no known high velocity gas. Nelson et 
al. (2000) claim that the high velocity components gen- 
erally account for a small fraction of the total flux in the 
[OIII] emission lines. In several cases they find clouds with 
multiple velocity components. 

4.1. Oxygen Lines: [OIII], [Oil], and [01] 

The rcsuhs for [O III]/H/3 and [O II] /W arc shown in 
Fig. 2. Different models are selected to fit the data as 
expected from the multiple structure of the observed lines 
(Kaiser ct al. 2000). As shown in Fig. 2 a, the data arc well 
reproduced by composite RD models, with the intensity 
of the ionizing radiation Fh decreasing by more than two 
orders of magnitude from the center towards the outskirts 
of the nuclear region. 

Focusing on [0III]/H/3, RD clouds (Fig. 2 a) with 
velocities of 200-300 kms~^ contribute preferentially in 
the inner 1 arcsec region, while those with velocities of 
100 kms~^ appear at larger distances from the center, 
in agreement with Winge et al. (1999). The contribution 
to [0III]/H/3 line ratios from RD clouds with Vs~ 500 
kms~^ is small. 

Regarding the SD clouds (Fig. 2 b and c), the high 
velocity ones (Vs> 500 kms""'^) are responsible for only a 
few percent of the central region emission, while a larger 
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Fig. 5. The distribution of the electron temperature 
and of the ionic fractional abundance of oxygen down- 
stream for a shock-dominated model (Fh=0) with Vs=700 
kms~^and no=700 cm~'^. The shock front is on the left. 



contribution at 3" from the center comes from the low 
velocity clouds (Vs< 500 kms~^). 

The situation is markedly different regarding 
[OII]/H/? (Figs. 2 d, e, and f) as it appears domi- 
nated by SD models. The trend of [0II]/H/3 is nicely 
explained by SD models with Vg between 400 and 700 
kms""'^ (thick line). 

Recall that the contribution of high velocity clouds to 
the [O II]/H/3 line ratio is due to the diffuse radiation 
generated at the high temperature post-shock zone reach- 
ing the low ionization zone. In order to illustrate this, the 
distribution of the temperature as well as the fractional 
abundance; of oxygen ions downstream arc plotted in Fig. 
5 for an SD model with Vs= 700 kms^^. 

Notice, however, that explaining the [0II]/H/3 line ra- 
tio only with SD clouds may not be acceptable in a global 
picture for the NLR, since RD clouds must be present 
and contributing to other emission lines. Indeed, the re- 
duction by a factor of 3 is not only dictated by the fit 
of the [0I]/H/3 data (Fig. 3 b), but is consistent with a 
general scenario. 

Regarding the [01]/H/3 line ratio, notice that SD mod- 
els overpredict the observed [O I]/H/3 data by a factor of ~ 
3. This indicates that SD models with Vs~ 700 kms~^ in 
the central region must be taken with a lower weight in 
an eventual averaged model. The shift of the SD models 
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dictated by the [0I]/H/3 ratios is consistent with all the 
other lines. 

4.2. Sulphur Lines: [SIII] and [SII] 

To further investigate NGC 4151 line emission we address 

a problem raised by Kraemer et al. (2000), namely, the 
overprediction of the [SIII] 9069+9532 /H/3 line ratio by 
photoionization models. This ratio is shown in Fig. 3. 

Different clouds with different Vs coexist in the center. 
The larger values can be fitted by RD models with Vs> 
500 kms~^. Notice, however, that RD models with Vs= 
500 kms"-'^ and log Fh< 10.5 ovcrpredict the data in the 
region beyond 1" at both sides of the center. Therefore, 
for consistency, they were taken out from all figures. 

There is no contribution from SD clouds to the 
[SIIIJ/H/? ratio in the central region (Fig. 3 e, f). On the 
other hand, a strong contribution to the [SII] /H/3 line ra- 
tios comes from the SD clouds with rather high velocities 
(700 kms~^) in the central region, and a not negligeble 
contribution of SD clouds with Vs= 1000 kms~^ (Fig. 4 
c), as for the case of [OIIJ/H/?. Here too the signature of 
the diffuse radiation generated at the high temperature 
post-shock zone is seen. 

From the [SII]6717/6731 diagrams of Nelson et al 
(2000, Fig. 10), the gas is rather clumpy and also indicates 
a decrease in density with distance. The data in the outer 
region (> 2") arc well explained by RD models with Vs= 
100 kms-i and uq = 100 cm'^, leading to [SII] 6716/[SII] 
6730 > 1. Density downstream decreases for lower shock 
velocities and lower preshock densities (see Contini & Al- 
drovandi 1986), therefore, [SII] 6716/ [SII] 6730 > 1 is ex- 
pected in the regions farther from the center. In fact, we 
have found that velocities generally decrease with distance 
from the center. However, some ratios of about 0.5 are 
given by Nelson et al. They are supported by the fact 
that some models with Vs of about 500 kms~"'^also fit the 
data beyond 1" (Figs. 4b and 4c). 

4.3. The [FeVII] line 

Fe coronal lines are usually underpredicted by pure 

photoionization models. In the case of NGC 4151, 
[FeVII] /H/? is underpredicted by a factor of ~ 3 by Krae- 
mer et al. (2000) and by Alexander et al. (1999). To ex- 
plain the discrepancy, an overabundance of Fe is often sug- 
gested, which is somewhat surprising since a fraction of Fe 
may be locked in grains, although sputtering is strong for 
small grains and high velocity shocks. 

[FeVII]/H/? values for RA. = 221° are shown in Fig. 
4 d, e, and f. The available data are scarce, hence fur- 
ther constraints to the models arc really not possible. The 
models used to fit the other lines are compared with the 
available data. As seen in Fig. 4 f high velocity models 
overpredict the [Fe VII] /H/3 line ratio. A better fit is ob- 
tained reducing the weight of the high velocity models 



by a factor of 2. This reduction does not change our dis- 
cussion above concerning the other lines, because models 
with Vs=1400 km s~^ contribute by no more than 20 % 
to [Oil] /H/3 and [SII] /H/3 line ratios. As we will sec over 
the next sections, a low contribution from high velocity 
clouds is also needed to explain the IR lines (§6) and the 
continuum SED (§7). Notice that the [Nc V] 3426 lines, 
also coming from the high ionization zone in the clouds, 
are overestimated by high velocity models (see §5). 

5. The Ne and Si infrared lines 

The ISO coronal line data reported by Sturm et al. (1999) 
are used to further constrain the models. Since the ISO 
aperture are of the order of several arcminutes, only the 
coronal lines are used. Lines from lower ionization stages 
may include the contribution from star forming regions in 
the galaxy. Nevertheless, they are included in the discus- 
sion in order to complete the data set with information 
from different ionization stages. 

Two Ne coronal lines are observed : [NeVI] 7.6 and 
[NeV] 14.3 in addition to two other strong lines, [Nelll] 
15.5, and [Nell] 12.8. Regarding Si, [Si IX] 3.94, [SiVII] 
2.48, and [Sill] 34.8 are detected. 

From all the models plotted in Figs. 2 to 4, we select 
those which could better reproduce the IR line ratios. The 
corresponding input parameters are listed in Table 1, as 
well as the weights (W(Ne)) adopted for the best average 
model (AV). The weights are relative to that of model 2. 

The weights adopted for the Si lines (W(Si)) are, how- 
ever, different because they account for Si depletion from 
the gas phase (and included in grains), which is different 
in each model. To account for the Si coronal line emission, 
Si/H has to be depleted by a factor of about 15. Moreover, 
the different depletion of Si in the different clouds may be 
indicating that the dust is not homogeneously distributed 
in the emitting clouds (see Kraemer et al. 2000). 

Notice that the weights are very different for different 
models because they must compensate the difference by 
many orders of magnitude of the line intensity fluxes, in 
order to obtain similar contribution to the emission lines 
from different models (CVOl). SD and RD clouds corre- 
spond respectively to models 2, 3, 4, 6, and models 1 and 
5. 

The observed IR emission-line flux as well as the ratios 

of the calculated to observed lines are listed in Table 2. 
This table is organized as follows: (a) Line fluxes (row 
1); (b) the results corresponding to the best fit model of 
Alexander et al. (1999) (row 2); (c) single-cloud results 
(rows 3 to 8); (d) the weighted average model AV (row 
9); (e) the contribution to the IR emission lines from each 
single-cloud model (rows 10 to 15). 

Actually, the line and continuum spectra must be mod- 
eled consistently, so, the results of the AV model in Table 
2 are cross-checked by the results of the continuum SED 
in Sect. 7, until the best tuning for both is achieved. 
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Table 1 

The models for the infrared lines 



model 1 2 3 4 5 6 



Because the theoretical fluxes are calculated at the 
nebula and the observed fluxes are measured at Earth, 
the ratio of the square distance of the galaxy to Earth 
(d^) to the square distance of the nebula to the galaxy 
center (R^) is given in the last column of Table 2. 

From Table 2, one sees that model 2 (SD, Vs= 300 
kms"'^) strongly contributes to all the lines. Gas is heated 
to temperatures between 1.3 10^ K and 3.7 10^ K for SD 
clouds with Vs of 300 kms~^ and 500 kms~^, respectively. 
Radiation dominated model 1 contributes mainly to [NeV] 
and [Sill] and model 5 to [Nelll] and [SiVII]. Model 5 is 
associated with a very high ionizing flux: 5 10^^ (cm~^ 
eV~^ at the Lyman limit). 

The high velocity clouds (model 6, Vs= 1500 kms~^) 
do have a low contribution. This contribution must remain 
low as to prevent the increase of [SilX] flux far beyond its 
observational value. 



dVR' 



10^ 

lO'' 

10^ 

10^3 

10« 



6. The UV line ratios 

Ultraviolet lines may provide important information 

about their origin in a gas cither heated by shocks or ion- 
ized by a strong radiation. The most significant lines are 
NV 1240, CIV 1550, and Hell 1640. Other strong UV lines 
are generally blended (OIV, SilV 1402, GUI], Silll] 1909, 
etc), and have not been used in the modeling. Because the 
reddening correction in the UV is important, we use the 
NV/GIV and Hell/CIV line ratios instead of intensities 
relative to H/3. 

The UV data correspond to P.A.= 221° observations. 

Model results are presented in Fig. 6 where South- West 
(SW) data (filled triangles) are separated from those form 
the North-East (NE) region (open triangles). 

The theoretical results indicate that the NE UV data 

can be explained by low velocity RD clouds (< 300 
kms~^) reached by an ionizing flux Fh< 10^^ units, while 



Vs(kms-i) 


100 


300 


500 


500 


500 


1500 


no (cm~^) 


100 


300 


300 


300 


300 


300 


log Fh 


11 








12.7 




D (10^^ cm) 


100 


100 


10 


100 


100 


1 


W(Ne) 


5.(-4) 


1. 


6.5(-4) 


7.9(-4) 


2. (-6) 


2.9(-8) 


W(Si) 


0.114 


1. 


1.4(-3) 


1.43(-3) 


1.4(-4) 


1.6(-6) 



Table 2 

The infrared lines relative to the observation data 



fluxes (obs)^ 


[NeVI] 


[NeV] 


[Nelll] 


[Noll] 


[SilX] 


[SiVII] 


[Sill] 


7.90 


5.50 


20.7 


11.8 


0.41 


1.2 


15.6 


model (calc/obs)^ 
model 1 (calc^/obs) 


0.5 


1.8 


0.7 


0.5 


0.9 


0.4 


3. 


2.30 


10.3 


2.39 


0.017 


0.064 


0.137 


80.4 


model 2 (calc/obs) 


1.20 


0.40 


0.50 


0.90 


17.6 


3.90 


6.6 


model 3 (calc/obs) 


1.80 


0.70 


6.80 


11.2 


205. 


11.1 


36.7 


model 4 (calc/obs) 


1.80 


0.70 


6.70 


11.0 


199. 


10.6 


27.7 


model 5 (calc/obs) 


1.98 


15.0 


27.5 


0.92 


3.70 


5.42 


0.11 


model 6 (calc/obs) 
model AY (calc'*/obs) 


2.05 


1.53 


0.16 


0.01 


1.04e5. 


240. 


0.0 


1.30 


1.24 


1.19 


0.92 


1.44 


0.80 


1.0 


model 1 % 


8.70 


42.8 


10.3 


0.09 


3.04 


11.9 


57.3 


model 2 % 


88.3 


32.3 


42.0 


97.7 


72.7 


29.4 


41.2 


model 3 % 


0.01 


0.0 


0.04 


0.08 


1.21 


0.12 


0.33 


model 4 % 


0.01 


0.0 


0.04 


0.09 


1.18 


0.11 


0.25 


model 5 % 


3.00 


24.9 


47.5 


2.05 


21.8 


58.5 


0.98 


model 6 % 


0.0 


0.0 


0.00 


0.0 


0.06 


0.0 


0.0 



^ in 10"" ergcm"^s"^ 

^ Alexander et al. 

^ calculated at the nebula 

* To calculate the averaged model from the models given in CVOl : 
Ei [(F(A)/F(H/3))cvoi F(H/3)ovoi)i W(model i)]/ F{X)obs 
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Fig. 6. He II 1640/C IV 1550 versus N V 1240/C IV 
1550. SW UV data and NE UV data correspond, repec- 
tively, filled and empty triangles. The RD results for Vs = 
100 kms-^and no=100 cm-'^ (dotted line), and Vs=200 
kms^^and no— 200 cm~'^ (dash-dotted line) are labelled 
by the values of log(FH), and the SD results (sohd line) 
by Vs in kms"^. 



the SW UV data come from the RD clouds reached by a 
stronger radiation field, with some contribution from SD 
clouds with Vs> 500 kms~^. However, these results are 
only indicative (see the discussion by Kraemer et al. 2000). 
It is well known that the Hell 1640 line is strongly depen- 
dent on the spectral index of the ionizing radiation, and 
the NV/CIV ratio depends on the adopted N/C abun- 
dance ratio. Regarding the Hell lines, if auv < 1-5 the 
theoretical results may vertically shift to the upper part 
of the diagram, and the derived velocities would be larger. 
Regarding the N/C abundance ratio, the cosmic values is 
adopted. 

For consistency, we present in Table 3 the ob- 
served/calculated UV-optical line ratios obtained with the 
same AV model as for the IR lines (see Table 2). To avoid 
the uncertainty of the dilution factor the ratios of the AV 
model are calculated between the line ratios to H/3. 

The fit is within 1.6 and the trend is always to overes- 
timate the values. This is consistent with reduction of the 
SD models which was adopted in Figs. 2-4 b, c, and e, f . 



7. The continuum 

To further constrain our model for NCC 4151, the contin- 
uum SED is analysed and modeled using the composite 
model derived from the Ne and Si lines in §5. 

Observations 

Gamma-rays to the optical continuum emission is 
taken from the 1993 December multiwavelength monitor- 
ing campaign of NGC 4151 (Edelson et al. 1996). The 
source was near its peak historical brightness during this 
campaign and showed the strongest variations of its con- 
timmm emission at medium energy X-rays (^1.5 keV) 
with amplitude variations of 24%, weaker variations (6%) 
at the gamma-ray energies, decreasing variation from ul- 
traviolet (9%) to optical (1%) and finally, not significant 
variation at the soft X-rays (0.1-1 keV). 

Assuming that the continuum emission variations are 
within a few percent level beyond optical wavelengths, the 
1993 data are combined with IR data taken at different 
epochs. Since the optical aperture used by Edelson et al. 
(1996) is 12x3 arcsec, near -IR fluxes from an aperture 
size of 10 arcsec is adopted when possible. Sources are as 
follows: 

The data for the near-IR come from the NASA Extra- 
galactic Data (NED), J, H and K data are taken from 
Balzano & Weedman (1981) and L and M data from 
Mc Alary et al (1979). The 10 microns emission is from 
Lebofsky & Rieke (1979) and corresponds to an aperture 
size of 6 arcsecs. In addition, to trace the non-stellar con- 
tribution at the near-IR, Kotilanien's et al. (1992) data 
within a 3 arcsec aperture are also considered for compar- 
ative purposes. 

For the far-IR region between 16 and 200 microns, 
ISOPHOT data from Perez-Garcia et al. (1998) are taken. 
Due to the large aperture size used in ISOPHOT, those 
data are integrated over the complete galaxy. 

For the radio data, values at 1.4 GHz, 4.85 GHz, and 
408 MHz were at the NED, and come from Becker et al. 
(1995), Becker et al. (1991), and Ficarra, Grueff & Tom- 
masctti (1985), respectively, while data at 8.4 and 5 GHz 
arc taken from Pedlar et al. (1993). 

Compared with other Seyfert galaxies, NGC 4151 is 
relatively weak in X-rays with L(2-10 keV) ~ 710^^ erg 
s^^ (Weaver ct al. 1994). The soft X-ray emission is ex- 
tended (Morse et al. 1995). Recent Chandra data resolve 
up to ^ 70 % of the 0.4-2.5 keV emission (Ogle et al. 
2000). This emission appears to be associated with the 
optical narrow line gas (NLR) extending asymmetrically 
to the South- West of the nucleus. 

On the nature of the soft X-ray emission 

Both, Weaver et al. (1994) and Ogle et al.(2000) pro- 
vide similar estimates for the plasma pressure of the hot 

gas of about 3 x cmr^ K . This is one order of magni- 
tude larger than the pressure derived from the cold NLR 
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Table 3 
The UV-optical lines 



[NcIV]2423 


[No V] 3426+ 


[NeIII]3869+ 


Hell 4686 


[OIII]4363 


[NII]6584+ 


1.77 


1.56 


1.81 


0.31 


0.35 


2.16 


0.4 


1.5 


0.8 


1.1 


0.7 


0.8 






0.5 


1.1 




1.8 


1.3 


1.5 


1.5 


1.1 


1.6 


1.4 



fluxes (obs)^ 
model (calc/obs)^ 

model (calc/obs)'' 
model AV (calc/obs) 



^ in 10"^^ ergcnr^s"^(Kraemer ct al 2000, Table 1, 0".l-0".3) 
^ Kraemer et al. (2000, Table 1, 0".l-0".3) 
^ Alexander et al. (1999) 

gas (Penston et al. 1990). The difference strongly argues 
against the soft X-ray emission being associated with the 
NLR confining medium (Weaver et al. 1994, Ogle et al. 
2000). 

Weaver et al. (1994) also found that the Fe K edge en- 
ergy in their 0.4 - 11 keV data (from the BBXRT mission) 
is inconsistent with an origin in a gas with the same ion- 
ization parameter as the low-energy absorber. This allows 
the authors to rule out a line-of-sight ionized absorber as 
the sole source of the soft X-ray excess in NGC 4151 (see 
also Contini & Viegas 1999 for the case of NGC 4051). 

George et al. (1998) propose a model for the 0.2-10 

keV spectrum where the underlying power law nuclear 
component is partially absorbed by an ionized absorber 
and partially scattered. However, it requires an additional 

component due to Brcmsstrahlung thermal emission from 
an extended photoionized gas at T'--^ 6 lO^K. 

Recent Chandra data provide direct evidence of X-ray 
line emission gas at T ~ lO'' K. The strength and ioniza- 
tion potential of the X-ray narrow emission lines indicate 
a composite spectrum in which both photoionization and 
shocks are at work (Ogle et al. 2000). The fair spatial as- 
sociation between the central NLR optical region and the 
soft X-ray emission indicated that both mechanisms are 
contributing to the multi-wavelength SED and to the line 
spectrum of NGC 4151 (cfr. Komossa 2001). Taking into 
account these observational constraints, the origin of this 
soft X-ray excess is evaluated in the next section assuming 
high velocity models dominated by shocks. 

Model Results 

NGC 4151 shows the characteristic of Seyfert galaxies 

rather than starbursts. Therefore, we consider that the 
continuum SED is Brcmsstrahlung radiation from clouds 
ionized and heated by the radiation from the AC and by 
shocks. We refer to previous studies (e.g Contini & Viegas- 
Aldrovandi 1990, Contini & Viegas 1991, Viegas & Contini 
1994, Contini & Viegas 2000, etc) on this subject which 
could lead to a better understanding of the results. Partic- 
ularly, dust and gas are coupled entering the shock front 
and mutually heat each other. The grains are coUisionally 
heated to the highest temperatures (> 300 K) leading to 
emission in the near-IR. On the other hand, heating by 



the central radiation is not eSicient enough, so dust in 
the radiation dominated zone does not reach such high 
temperatures. 

Before discussing model results, let us notice that be- 
low the Lyman limit, up to say 0.2 keV, we have very little 
observational information: this is the "unknown window" 
where observations arc difficult because of the heavy ab- 
sorption by our Galaxy. However, in the 0.2 to 2 keV re- 
gion, NGC 4151 shows extended emission, spatially corre- 
lated with the optical NLR gas. If due to Brcmsstrahlung, 
the tail of this X-ray emission should somehow show up 
in the far UV. However, no detected extended emission 
appears in the UV HST images of NGC 4151 (Boksen- 
berg et al. 1995). Therefore, the emission from the NLR 
gas should be lower than that of the UV nucleus of NGC 
4151. We will use this observational fact as a constraint 
in the proposed modeling in the sense that the observed 
nuclear UV emission can only be associated with model 
emission within a nuclear region of outmost 3 pc size. 

Observational data and model results are plotted in 
Figs. 7 and 8, which refer to models and relative weights 
adopted to fit the Ne infrared lines (Table 1, row 5). Notice 
that Brcmsstrahlung peaks at high frequencies depending 
on Vs(see Viegas & Contini 1994). 

The large geometrical thickness of the emitting clouds 

and the high densities downstream due to compression 
(n/uo > 10, depending on Vg, uq, and Bq) lead to high 
optical thickness of the emitting gas with column densities 
of the order of 10^^-10^^ cm^^. This is within the range of 
column densities estimated from the X-ray data, ~ 5 10^^ 
cm~^ (Weaver et al. 1994, George et al. 1998, and Ogle et 
al. 2000). Accordingly, the emitted radiation gets absorbed 
in the cloud itself between ~ 13.6 eV and 500 eV, hence 
the emission gap in this region shows up in Figs. 7 and 
8. The high velocity model 6 has however a lower column 
density (< 10^^ cm~^) because is radiation-bounded. In 
fact, due to strong compression, the temperature rapidly 
decreases downstream at a distance < 10"'^^ cm from the 
shock front. Thus, the contribution of model 6 to the SED 
extends into the far UV range (Fig. 7 long-dashed line). 
However, the X-ray data in the frequency range 10"'^^ and 
10^^ Hz impose a limit on the relative constribution of this 
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Fig. 7. Different components contributing to the SED of 
the continuum of NGC 4151, corresponding to the results 

of several models, scaled following the weights used to fit 
the IR lines. The models are the following: 1 (dash-dotted 
line), 2 (dotted line), 3 (short dash-dotted line), 4 (short- 
dashed line), 5 (solid line), 6 (long dashed line). The flux 
from the AC in the far X-rays, which is seen through 
the clouds, is represented by a thick dotted line, while 
the continuum from the old stellar population is indicated 
by a thick long dash-dotted line. The theoretical results 
are compared to observational data taken from the NED 
(filled squares) and from other sources described in Sect. 
7 (open squares). 



model 6 (see below). All together, the fit of the continuum 
at the UV range is obtained with model 2. The dust-to- 
gas ratio that characterizes this model is d/g = 10^^"'. 
The other models (models 3, 4, 5, and 6) underpredict the 
data. 

The summed SED which better explains the obser- 
vations is shown in Fig. 8. The main components are 
summed separately: radio emission, reradiation by dust in 
the IR, gas Bremsstrahlung, and the flux from the active 
source. Moreover, as was found for many Scyfcrt galaxies 
(see Contini & Viegas 2000) a black-body with tempera- 
ture of ^ 3000 K is used to roughly represent the optical 
continuum due to the stellar population (thick solid line). 

Notice that we are interested in the bulk contribution 
of the stars rather than the detailed shape of their contin- 
uum emission. Thus, although we could use the spectrum 
of an early-type galaxy or the spectrum of the central re- 



Fig. 8. The continuum spectrum of NGC 4151. The thin 
solid lines correspond to the summed theoretical contribu- 
tions from dust emission, gas Bremsstrahlung and radio 
emission as shown in Fig. 7. The thick solid line corre- 
sponds to the summed contributions from dust emission, 
gas Bremsstrahlung and the old stellar population contin- 
uum, while the dotted line corresponds to the AC con- 
tinuum seen through the clouds. The observational data 
are taken from the NED (filled squares) and from other 
sources described in Sect. 7 (open squares). 

gion of a Sa galaxy to get a more realistic fit of the stellar 
contribution, it would not change our conclusions. 

For instance, Malkan & Filippenko (1983) find that 
about 1/10 of the emission at 5500 land ~ 1/3 at 9000 
Aare due to stars (in aperture of 10" , just the same aper- 
ture of the optical and IR data used in this paper), namely 
they estimated the stellar emission to be 14 mJy at 5500 
A. We have added the corresponding value to Fig. 7 (star) . 
Accounting for the errors in the observations and models, 
notice the good agreement (within 30%) with the black 
body emission curve. 

The bump in the IR is rather wide, representing the 
sum of contributions from different models. In previous 
papers (e.g. Contini et al. 1999b, Appendix A) it was ex- 
plained that the bump in the IR due to each model de- 
pends on the shock velocity. Particularly, the high velocity 
model (model 6) has its maximum at the near-IR. To fit 
the near-IR data, a high d/g (2 10~^^) is assumed in model 
6, indicating that high velocity material is rather dusty. 
A shock velocity of 1500 kms"-'^ produces a post shock 
region at a temperature of 3.4 10^ K, in agreement with 
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the hot plasma temperatures derived from the X-ray data, 
indicating that high velocity material is rather dusty. 

Shocked clouds with Vs~ 700 kms~^ produce lower 
temperatures, ^ 10^ K, which are still within the range 
of values estimated from X-ray data (e.g. George et al. 
1998). These velocities are also revealed by optical lines 
(§4). 

Yet, the contribution of a model with Vs= 1500 
kms~^ in the SED should be taken with a relative small 
weight (< 2.9 10^^), because the data in the soft X-ray 
range constrain the model. As a consequence, the contri- 
bution of this model to the IR lines is very low. Thus, a 
signature of these high velocity clouds is not expected to 
be seen in the infrared emission-line spectrum. 

In the radio range. Pedlar et al. (1993) derive an aver- 
age spectral index for the entire source of a = -0.87 from 
their 8.4 and 5 GHz data. They interpreted the emission 
as due to thermal free-free emission from the NLR. In the 
present modeling, the radio data are instead well explained 
by synchrotron radiation with spectral index -0.75, gener- 
ated by Fermi mechanism at the shock front (Bell 1978). 
Notice that self- absorption of the Bremsstrahlung emis- 
sion at these radio wavelengths may be important, hence 
its contribution is negligeble (see, for instance Contini et 
al 1998b). 

8. Discussion and Conclusions 

In this paper we have modeled the narrow emission-line 
and continuum spectra of NGC 4151 with particular at- 
tention to the large range of velocities indicated by the 
line profiles. 

The analysis of the line profiles is complex and clouds 
in many different conditions can contribute to each line. 
By modeling the NLR to analyse the high spatial resolu- 
tion data (Nelson et al. 2000, Kraemer et al. 2000), it is 
found that the contribution of each cloud to a given line 
may show a large variation from line to line. 

In the central region high velocity clouds revealed by 
the observations (Hutchings et al. 1999, Winge et al. 1997, 
etc) are shock dominated (SD). There is a strong contribu- 
tion of SD clouds with shock velocity Vs~ 700 kms~^ to 
[Oil] 3727/H/3, and [01] 6300/H^. Radiation dominated 
(RD) clouds with Vg^ 500 kms"'^ are necessary to explain 
[SIII]/H/? and [OI]/H/?, while clouds with low (100-200 
kms~^) and reached by a relatively strong ionizing radia- 
tion (log Fh'- 11) contribute to the [0III]/H/3, [0II]/H/3, 
and [SIIJ/H/3 line ratios. 

Beyond 2" from the center, RD clouds, photoionized 
by a weaker radiation flux are responsible for the observed 
line emission. 

Quantitative modeling of the Ne and Si infrared lines 

from different ionization levels leads to more results. No- 
tice, however, that the observations give the integrated 



flux from all the galaxy, so that it is not possible to un- 
derstand which type of clouds prevail in different regions. 

SD clouds with Vs= 300 kms~^strongly contribute 
to all Ne lines corresponding to different ionization lev- 
els (Table 2), while RD clouds with velocities of 500 
kms~^ and reached by a strong radiation flux (log Fh= 
12.7) contribute particularly to the [Nelll] and [SiVII] IR 
lines. Due to the high postshock temperature, high veloc- 
ity clouds (Vs= 1500 kms~^) would overpredict the [SilX] 
3.94 line. Therefore, the contribution of these clouds is low. 
A large contribution to [NeV] and [Sill] lines come from 
low velocity (100 kms~^) clouds reached by an ionizing 
radiation characterized by log Fh= 11- 

Modeling of the IR lines also shows that silicon is de- 
pleted by a factor of 15 because included in dust grains, 
while N/C abundance ratio is found compatible with cos- 
mic values from the modeling of the UV lines. 

Comparing the results of the present work with those 
of other authors, notice that our results are slightly better. 
Indeed, better fits could be obtained with a different choice 
of weights. However, the results presented in Tables 2 and 
3 are consistent with the fit of the continuum SED (§7) 
and were chosen accordingly. More particularly, it seems 
that the three Ne emission lines are higher than observed 
by factors of 1.3 - 1.5. This may be an indication that 
the Ne abundance we have used is a factor of 1.3-1.5 too 
high. So, decreasing the Ne/H abundance by 1.3, we get a 
better agreement. The Hell 4686 line is slightly higher, but 
this is mainly due to the power law above 54.4 eV, which 
we (and the others) could have taken somehow too flat. 
The only problem is the [0III]4363 line. Because oxygen 
is a coolant, changing the O abundance may not solve the 
problem. However, a factor of 1.6 higher is still within the 
limit of a factor of 2 proposed by Alexander et al., who do 
not give their results for [OIII] 4363. 

Notice that we based our models on IR lines and we 
compared with Alexander et al. (1999). So, we should re- 
fer only to them. The results of Kraemer et al in row 2 
of Table 3 are better than ours only for [OIII] 4363 and 
[Nelll]. Actually, they come from observations in the re- 
gion between 0.1" and 0.3" SW, whereas we average on 
all the regions. 

The analysis of the continuum SED leads to the fol- 
lowing results: 

The high velocity material is very dusty (d/g~ 10~^^). 
The radio emission is synchrotron created by Fermi mech- 
anism at the shock front. Moreover, shocks are important 
to explain the soft X-rays: shocks with velocities of 1500 
kms~^ produce a post shock region at a temperature of 
3.4 lO'^ K. This value is in agreement with the temperature 
found by Weaver et al. (1994) and Ogle et al. (2000) for 
a non-equilibrium plasma from X-ray observations, while 
clouds with Vs~ 750 kms~^ produce temperatures of 8. 
10^ K (George et al. 1998). 
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An estimate of several physical quantities applying to 
the central AGN region can be derived from the proposed 
modeling. 

If d is the distance from Earth (19.8 Mpc) and d^/R^ 
= 10^ (see Table 2), R ~ 0.626 kpc is the average distance 
of the emitting nebula to the AC. Adopting an average 
downstream density n=10''' cm^'', and D=10"'^* — 10^^ cm, 
the calculated mass is M ~ 3.5 10^~^ ff Mq , where ff is 
the filling factor which is likely to be less than unity. In 
this case our estimation of the NLR mass is lower than 
the values quoted by Ulrich (2000) for an emitting region 
closer to the center, i.e. 10^ Mq within 40 pc and 5 10^ Mq 
within 12 pc (0.15"). The corresponding calculated central 
source luminosity is about 1.6 10'''^ erg s~^. The average 
kinetic energy of the NLR is about 3 10^^ ff ergs, assuming 
an averag e velocity of the NLR clouds of 300 km s" ^ . This 
would imply that the time scale for the AGN phase in 
NGC 4151 should be larger than 5 10^ years. 

Summarizing, we provide a self-consistent modeling of 
the multiwavelength line spectrum and SED of the nuclear 
region of NGC 4151. 

The modeling is based on the coupled effect of shocks 
and photoionization operating in the narrow emission line 
gas. As such, it implies the existence of emitting clouds 
with velocities and densities in a large range. An emitting- 
cloud distribution in the velocity/density space is precisely 
what is revealed by the HST observations of this galaxy. 
This together with the fair fitting of the line and contin- 
uum spectra obtained in this work reinforces our hypothe- 
sis that shocks and photoionization are effectively coupled 
in the NLR of AGN. 
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